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We report the results of the first search for gravitational waves from compact binary coalescence using data 
from the Laser Interferometer Gravitational-wave Observatory (LIGO) and Virgo detectors. Five months of 
data were collected during the concurrent S5 (LIGO) and VSR1 (Virgo) science runs. The search focused on 
signals from binary mergers with a total mass between 2 and 35 M©. No gravitational waves are identified. 
The cumulative 90%-confidence upper limits on the rate of compact binary coalescence are calculated for non- 
spinning binary neutron stars, black hole-neutron star systems, and binary black holes to be 8.7 x 1CP 3 yr -1 , 
2.2 x 1CT 3 yr“ 1 Li 0 1 , and 4.4 x 10 -4 yr _1 L^ 0 1 respectively, where Lio is 10 10 times the blue solar luminosity. 
These upper limits are compared with astrophysical expectations. 

PACS numbers: 95.85.Sz, 04.80.Nn, 07.05.Kf, 97.80.-d 


I. INTRODUCTION 

The coalescence of a stellar mass compact binary is ex- 
pected to produce gravitational waves detectable by ground- 
based interferometers. Binary neutron stars (BNS), binary 
black holes (BBH) and black hole-neutron star binaries 
(BHNS) can spiral together to produce signals in the fre- 
quency band where the Laser Interferometer Gravitational- 
wave Observatory (LIGO) [1] and Virgo [2] detectors are most 
sensitive (40-1000 Hz). 

LIGO was collecting data at the Hanford, Washington and 
Livingston, Louisiana sites as part of its fifth science run (S5) 
(4 November 4 2005 - 30 September 2007) when the first sci- 
ence run (VSR1) began at the Virgo detector in Cascina, Italy 
on 18 May 2007. During VSR1 the Virgo detector operated 
at reduced sensitivity since its commissioning was still incom- 
plete. LIGO data collected before 18 May 2007 were analyzed 
separately and upper limits on the rate of gravitational waves 
from binary inspirals were reported in Refs, [3-5]. 

Here we describe the results of the first joint search for 
gravitational waves from compact binary coalescence with 
LIGO and Virgo data. This search covers gravitational waves 
from binaries with a total mass between 2 M 0 and 35 Mq and 
a minimum component mass of 1 Mq. This analysis is based 
on the same methods as the S5 LIGO-only searches [4, 5]. 
Since the analysis is considered integral in preparing for fu- 


ture joint searches in LIGO and Virgo data, further develop- 
ments were performed to integrate Virgo into the pipeline, 
even though VSR1 data had limited sensitivity when com- 
pared with LIGO’s S5 data. No gravitational-wave signal is 
identified and upper limits are calculated. 

In section n, we describe the data used in this analysis. The 
data reduction pipeline is explained in section m and ends 
with a description of the detection statistic. The results and 
upper limits appear in sections IV and V. Details of a self- 
imposed blind injection challenge are given in Appendix A. 


n. DATA QUALITY 

The detectors are referred to as HI (Hanford 4 km), H2 
(Hanford 2 km), LI (Livingston 4 km), and VI (Virgo 3 km). 
Data from LIGO and Virgo are recorded in the same format, 
making it easier to run the LIGO pipeline on the additional 
detector. The relative sensitivities of these detectors can be 
assessed with horizon distance, the distance at which an opti- 
mally located, optimally oriented binary would produce trig- 
gers with a signal-to-noise-ratio (SNR) of 8 in the detector. 
When averaged over the duration of the search, the horizon 
distances for a 1.4, 1.4 Mq BNS system are approximately 
37, 16, 32, and 8 Mpc for HI, H2, LI, and VI respectively. 
See Figure 1 for the horizon distance in each interferometer 
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as a function of the total mass of the binary system. 


HI. THE DATA ANALYSIS PIPELINE 



FIG. 1 : The average inspiral horizon distance over the run is shown 
as a function of the tot i mass of the binary system for each inter- 
ferometer. The error bars indicate variation over the duration of the 
run. 


The detectors are very sensitive to their environments and 
fall out of science mode when disturbed, meaning that they 
are temporarily not recording science-quality data. Because 
the data streams from each detector are not continuous, differ- 
ent combinations of detectors may be taking data at any given 
time. As we describe in section HI A, we require time co- 
incidence to identify possible gravitational waves and hence 
we only analyze the data when at least two detectors are 
operating. There are eleven combinations for what we de- 
fine as analysis time: H1H2, H1L1, H1V1, H2L1, H2V1, 
L1V1, H1H2L1, H1H2V1, H1L1V1, H2L1V1, H1H2L1V1. 
Analysis time indicates that the listed detectors are collecting 
science-quality data. Because HI and H2 are co-located, cor- 
related noise leads to poor background estimates and hence 
H1H2 time was rejected. 

A number of quality criteria were established before and 
during the run to reject times when the data are unreliable, 
either due to instrumental problems or external factors. See 
Appendix A of Ref. [4] for a more thorough description of 
the veto categories we use. We do not analyze data rejected 
by Category 1 vetoes because it indicates severe problems. 
Category 2 vetoes remove artifacts with well understood ori- 
gin and coupling. Category 3 vetoes are based on statistical 
correlations, and Category 4 vetoes are least serious and only 
used in the candidate follow up procedure. We make our can- 
didate event list and perform the upper limit calculation with 
data that passes Category 1+2+3 vetoes. We also look for 
loud candidates with significantly low false alarm probabil- 
ity that occur during times rejected by Category 3 (but that 
pass Category 1+2 vetoes). When only Category 1+2 vetoes 
are applied, 115.2 days of data are analyzed; when Category 
1+2+3 vetoes are applied, 101.1 days of data are analyzed. 
Our two most sensitive detectors (HI and LI) were simulta- 
neously running during 68% of this time. 


The data processing is performed in a similar manner to the 
S5 LIGO-only analyses [4, 5], although the addition of Virgo 
to the pipeline led to enhancements in the ranking method for 
candidates. Due to long-term variations in detector perfor- 
mance, data are analyzed in one-month blocks of time in order 
to obtain more accurate background estimates. There are four 
approximately 30-day blocks and one 19-day block and each 
time period is analyzed with an identical pipeline. The results 
of these five periods are combined with previous analyses into 
one set of upper limit statements. 


A. Overview of Pipeline 


As described in the LIGO-only searches [4], the analysis 
begins with four separate data streams, one from each detec- 
tor. We construct template banks [6] of non-spinning post- 
Newtonian waveforms [7-17]. These templates cover a range 
of binary mass combinations, (mi, m 2 ). The single-detector 
data are match filtered with the templates and the resulting 
triggers pass to the next pipeline stage if they exceed an SNR 
of 5.5 [18]. Because the background does not follow a Gaus- 
sian distribution, the false alarm rate (FAR) is quite high in 
single-detector data. To reject noise artifacts, we use signal- 
based vetoes [19, 20], including a x 2 test [21] and require 
triggers from different detectors to be coincident in time and 
mass parameters [22]. We define event type as the combina- 
tion of detectors contributing to a given coincident trigger. A 
double coincidence trigger can occur during double, triple, or 
quadruple analysis time, while a quadruple coincidence can 
only occur during quadruple analysis time. We apply consis- 
tency tests on the coincident triggers. For example, since HI is 
about twice as sensitive as H2, any coincidence that includes 
an H2 trigger, but not an HI trigger when HI was collect- 
ing data, is rejected. The remaining triggers are ranked based 
on an estimate of their likelihood of being a true signal or 
background. Any candidate events that stand out significantly 
above the background are followed up with a more detailed 
study of the triggers and detector conditions at the time of the 
event [23]. 

The background for the search is estimated by time-shifting 
the data from the different detectors. The time shifts are larger 
than the light-travel time between any pair of detectors, there- 
fore any observed coincidences in this data are accidental. The 
LI and VI data streams are shifted in increments of 5 and 
15 seconds, respectively, while the HI and H2 data streams 
are held fixed with respect to each other. This is because HI 
and H2 are co-located, and noise from environmental distur- 
bances is correlated in these interferometers. For this same 
reason, the background for H1H2 triggers can not be reliably 
estimated. H1H2 triggers are excluded from the calculation 
of the upper limit, but the loudest are followed up to ensure 
exceptional candidates are not missed. 
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B. Parameter Choices and Timing 

Many analysis parameters are determined at the onset of 
the analysis based on known properties of the individual de- 
tectors. LIGO data is analyzed above 40 Hz, and templates 
for the LIGO detectors cover a region with total masses be- 
tween 2 Mq and 35 M©. Virgo data quality information 
is best in the high frequency region, therefore the low fre- 
quency cutoff is set to 60 Hz for Virgo data. Consequently, 
the Virgo template bank is constructed to cover only the BNS 
mass region, with a minimum total mass of 2 M© and max- 
imum chirp mass of 2.612 Mq (where chirp mass is M c = 
((mim 2 ) 3 /(m 1 + m 2 )) 1 / 5 and ra x = m 2 = 3 Mq). 

When optimally tuned, veto cuts and consistency tests re- 
move a significant number of background triggers while hav- 
ing minimal effect on the detection efficiency for simulated 
signals. With the addition of a fourth detector, the tuning was 
revisited. In the process of tuning we set the appropriate pa- 
rameters for Virgo and verified that the corresponding param- 
eters for the LIGO detectors did not need to be changed from 
those used in S5 LIGO-only analyses. 


C. Detection Statistic 

In Refs. [4, 5], coincident triggers that survived all veto cuts 
and consistency tests [19] are ranked according to their com- 
bined effective SNR, p c , first used as detection statistic in the 
analysis of data from the S3 and S4 LIGO science runs [3]. 
The combined effective SNR statistic is based on the standard 
SNR, but it incorporates the value of the x 2 test into its def- 
inition [21]. Its effect is to assign a lower detection statistic 
to those coincident triggers that have high values of the x 2 > 
indicating that they are less consistent with the expected grav- 
itational waveform. Further details concerning construction 
of the combined effective SNR can be found in Appendix C 
of Ref. [4]. 

As observed in previous LIGO analyses, both the total 
rate of triggers and their distribution over effective SNR vary 
strongly with total mass. Variation also exists for each event 
type across different analysis times. Additionally, one should 
consider significant differences in detector sensitivities, for 
example the H2 and VI detectors are much less sensitive than 
either HI or LI. As a result, some analysis times are more ef- 
ficient in detecting gravitational waves than others. Within a 
specific analysis time, certain event types are more likely to be 
associated with a gravitational-wave event. Hence we specifi- 
cally distinguish all of the possible combinations of event type 
in analysis times. 

In order to account for variation in background rates and 
differences in the sensitivity of the detectors, we implemented 
the following post-processing algorithm. First, coincident 
triggers that survive the main pipeline are clustered such that 
only the trigger with the highest combined effective SNR 
within a 10 s window is kept. Then clustered triggers are sub- 
divided into categories by analysis time, event type, and mass. 
Based on regions of similar background behavior, we define 
three mass bins: 0.87 < M c /Mq < 3.48, 3.48 < M c < 7.4, 


and 7.4 < M c < 15.24. These correspond to equal mass bi- 
naries with total masses of 2 M© - 8 M©, 8 M© - 17 M© and 
17 Mq - 35 Mq. For every trigger in each category, using 
our estimate of the background (time-shifted data), we cal- 
culate the rate, i2 0 (p c ,m, /?), of background triggers with 

combined effective SNR greater than or equal to that of the 
trigger. The mass bin is indicated by m, while a and 0 are the 
event type and analysis time. Next, we introduce efficiency 
factors that estimate the probability of detecting a signal with 
a given combination of detectors in a specific analysis time. 
Virgo only has templates covering the BNS mass space, there- 
fore in the calculation of the efficiency factors we use a pop- 
ulation of simulated BNS gravitational-wave signals injected 
into the data. This procedure accounts for most of the effects 
introduced by variations in the detector sensitivities. Because 
the population of simulated signals is distributed uniformly in 
inverse distance, a reweighting is necessary. The efficiency 
factors are defined as: 


E D 3 . 

found mi /1 N 

e(a,p) = -= — -g-A (1) 

2^all 

The numerator is a sum of all injections found for that partic- 
ular a and 0. The denominator sums over all injected signals 
during a particular analysis time, 0. is the injected dis- 
tance to the binary. 

Finally we define the detection (or ranking) statistic, L, for 
the search to be 


L(pc,m,a,P) = In 


e(Q./3) 

R 0 (p c ,m,a,/3) 


( 2 ) 


For a gravitational- wave detection, a candidate is expected 
to have an L value significantly larger than the background. 
We have tested this algorithm on simulated signals and find 
that it results in substantial increase in overall efficiency of 
the search. 


IV. RESULTS 

A list of the loudest events is generated after Category 
1+2+3 vetoes are applied. However, in order not to unneces- 
sarily dismiss a possible detection, we also look for any loud 
candidates that might have occurred when a Category 3 veto 
was active (times that pass only Category 1+2 vetoes). Candi- 
dates from these times may still stand above the background, 
but must be closely studied to differentiate them from the el- 
evated background noise that the Category 3 veto is intended 
to remove. 


A. Results from Times that Pass Category 1+2+3 Vetoes 

After Category 1 +2+3 vetoes are applied, we find no events 
with a detection statistic significantly larger than the back- 
ground estimation. In Figure 2, the data are overlaid on the 
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background. The inverse false alarm rate is calculated with 
detection statistic L defined by Eq. (2). The loudest trigger in 
the five-month span is an H1L1 coincidence from H1H2L1V1 
time with a false alarm rate of 19 per year. As 0.28 yr was 
searched, this is consistent with the background expectation. 

However, as seen in Figure 2, there are fewer foreground 
triggers than the mean background. While the foreground lies 
within the 2 N 1 / 2 uncertainties, we performed a series of tests 
to exclude an error in the analysis or a bias in the way the 
foreground was handled with respect to the background. We 
ran the analysis on simulated Gaussian noise data. No deficit 
of foreground triggers in the tail of the distribution was found, 
which suggests that there is no problem in the analysis pro- 
cedure or codes. We studied how the data quality and ve- 
toes were applied in the analysis and found no error. We also 
changed the segmentation of the data and observed that the 
foreground events were shifted within the expectation for ran- 
dom fluctuations. Thus, we conclude that the results are con- 
sistent with a fluctuation of the foreground compared to the 
background. 



FIG. 2; A cumulative histogram of the inverse false alarm rate using 
L as the detection statistic. The data are represented by the triangles 
and each gray line represents a background trial made from time- 
shifting the data against itself. The darker and lighter shaded regions 
denote A 1/2 and 2N 1/2 errors, respectively. The data combine trig- 
gers from all five LIGO- Virgo months when Category 1+2+3 vetoes 
are applied. 


B. Results from limes that Do Not Pass Category 3 Vetoes 

An event list was generated for times when Category 3 ve- 
toes were active, meaning that the events were only able to 
survive the Category 1+2 vetoes. Only one event is incon- 
sistent with the estimated background. That sole significant 
candidate, an H1H2L1 triple coincidence, is a hardware injec- 
tion, part of a blind injection challenge. During four months 


of S5/VSR1, the LIGO and Virgo Collaborations agreed that 
simulated signals would be inserted into the LIGO-only data 
without the search groups knowing the time or number of in- 
jections and their parameters. This was an exercise to test the 
effectiveness of the search procedures and all blind injection 
triggers are removed from the results presented in this publi- 
cation. 

This sole candidate corresponds to a blind injection signal 
that was injected into the LIGO data during a time of high 
seismic activity at low frequencies at the LIGO Livingston 
Observatory. A Category 3 veto rejected this time period 
and hence this blind injection signal was not identified in the 
Category 1+2+3 event list. Unfortunately, the parameters of 
the blind injection challenge were revealed after the Category 
1+2+3 event list was produced, but before we looked for sig- 
nificant candidates that might have been removed by Category 
3 vetoes. Hence, the follow up procedure for significant can- 
didates was not exercised until after the injection parameter 
were known. Detailed investigations related to the blind in- 
jection challenge are described in Appendix A. 


C. Results for H1H2 Double Coincidences 

Although we do not have reliable estimates of the detection 
statistic for H1H2 events, we did look for interesting H1H2 
candidates and found one that passed Category 1+2+3 ve- 
toes. It corresponds to the same blind injection mentioned 
earlier. When the candidate was vetoed in LI, it became an 
H1H2 double candidate (see Appendix A). No other interest- 
ing H1H2 candidates are identified. 


V. UPPER LIMITS 

Other than the blind injection candidate, no significant can- 
didates are identified after Category 1+2+3 vetoes are applied 
or when Category 3 vetoes are disregarded. We calculate up- 
per limits on the rate of compact binary coalescence for the 
following astrophysical objects after Category 1+2+3 vetoes 
are applied: BNS [mi = m 2 = (1.35 ± 0.04) BHNS 
[mi = (5 ± 1) Mq, m 2 = (1.35 ± 0.04) M 0 ], and BBH 
[mi = m 2 = (5 ± 1) Mq]. We also present upper limits as 
a function of the total mass of the binary and as a function of 
the black hole mass for BHNS binaries. 

The upper limits are reported for both non-spinning and 
spinning objects in Table I. Only non-spinning templates are 
used in this search, so there is an additional loss of efficiency 
associated with spinning waveforms that leads to slightly less- 
constrained upper limits in the spinning case. The results are 
reported as a rate in units of number per Lio per year, where 
Lio is 10 10 times the blue solar luminosity, which is expected 
to be proportional to the binary coalescence rate [24]. The 
horizon distance listed in Table I is approximated for the HI 
or LI detector and is a good estimate of the sensitivity of the 
search. 

We calculate our upper limits using the loudest event from 
the search, as described in Ref. [25]. In this method, the pos- 
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tenor distribution for the rate depends on two quantities, Cl 
and A, that are functions of the loudness parameter, x. In our 
experiment, x is the inverse false alarm rate of the loudest ob- 
servec event according to the detection statistic in Equation 
2. A is a measure of the likelihood of detecting a single event 
with loudness parameter, x , versus such an event occurring in 
the experimental background. Cl is the cumulative luminosity 
of sources that produce signals that are louder than x . Assum- 
ing a uniform prior, the posterior distribution for the rate of 
coalescence is given by: 


p ( mIC L: T, A) = ^ (1 + pC L TA) e-^ T , (3) 

where p is the rate and T is the analyzed time. In general, A 
is given by [25]: 



where Pq is the probability of obtaining zero background 
events lcuder than x for the given search and observation time. 

The cumulative luminosity measures how many potential 
sources we can detect with this search, based on the blue 
light luminosity of galaxies. To find the cumulative luminos- 
ity, we take the product of the detection efficiency, calculated 
as a function of mass and distance, and the luminosity from 
galaxies in the catalog [24] and integrate over distance. We 
marginalize over our uncertainties when calculating the cu- 
mulative luminosity using the values given in Table I. These 
include detector calibration, Monte Carlo error, distances and 
luminosities given in the galaxy catalog and inaccuracies in 
the template waveforms [26]. The results from all five months 
and the prior S5 results [4, 5] are combined by taking the prod- 
uct of their posterior distributions calculated with uniform pri- 
ors as in Equation 3. Figure 3 shows the probability distribu- 
tion from the combined data for the rate of BNS coalescence. 

When spin is neglected and the priors from previous LIGO 
searches are used, the upper limits on the rate of compact bi- 
nary coalescence are 

7 ^ 90 %, bns = 8.7 x 10^ 3 yr^ 1 Lio _1 (5) 

7^90%, bhns = 2.2 x 10“ 3 yr l Lio 1 (6) 

^• 90 %,bbh = 4.4 x 10 _4 yr“ 1 Lio~ 1 , (7) 

which are consistent with upper limit estimates based solely 
on the sensitivity and observation time of the detectors [24]. 

Astrophysical observations of neutron stars indicate that 
their spins will be too small to have a significant effect 
on binary neutron stars (BNS) waveforms observable by 
LIGO [27, 28], hence we do not report upper limits for spin- 
ning BNS systems. However, we do consider spin effects on 
the upper limit for BHNS and BBH systems. The black hole 
spin, S', must be less than Gm 2 /c . We sample from a uni- 
form distribution of possible spin values in order to simulate 
the effect of spin on our ability to detect the binary system. 



Rate (yr 1 £, 10 1 ) 

FIG. 3: The posterior probability distribution for the rate of non- 
spinning BNS coalescence. The results of all previous LIGO 
searches are included in the plot as the prior, labeled as S5 LIGO- 
only Months 0-18. Each of the five LIGO- Virgo month results was 
combined with the prior to obtain the combined posterior, shown as 
the solid black line. 

With black hole spin included, the upper limits on the rate of 
compact binary coalescence are 

7 ^ 90 %, bhns = 2.7 x 10 3 yr 1 L 1 o 1 (8) 

7 ^ 90 %, BBH = 5.3 X 10 -4 yr -1 Lio -1 . (9) 

We also produce two sets of upper limits as a function of 
mass. The BBH upper limit shown in Figure 4 assumes a uni- 
form distribution of the component mass. The BHNS upper 
limit is shown as a function of black hole mass, assuming a 
fixed neutron star mass of 1.35 M 0 . 

VI. CONCLUSIONS 

We searched for gravitational waves from compact binary 
coalescence in the mass region 2 M 0 to 35 M 0 . Over 101 
days of coincident data were collected during the end of the 
LIGO S5 and Virgo VSR1 runs, making this the first joint 
search for gravitational waves from compact binaries with 
LIGO and Virgo data. The LIGO data analysis pipeline was 
augmented to handle the extra complexity of four detectors 
and a larger number of coincidence categories. Although no 
gravitational-wave candidates are identified, upper limits on 
rates of binary coalescence are established. The upper limits 
improve when combined with the previous LIGO-only results. 
These upper limits are still more than an order of magnitude 
larger than optimistic astrophysical expectations [29]. Hard- 
ware upgrades after S5 and VSR1 completed should yield bet- 
ter sensitivity in future searches. With the advent of three-site 
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BNS 

BHNS 

BBH 

Component Masses (M©) 

1.35/1.35 

5.0/1.35 

5.0/5.0 

Horizon Distance (Mpc) 

-30 

-50 

-90 

Cumulative Luminosity (lao) 

370 

1600 

8300 

Calibration Error 

13% 

14% 

14% 

Monte Carlo Error 

17% 

17% 

18% 

Waveform Enror 

19% 

18% 

16% 

Galaxy Distance Error 

-16% 

-13% 

-13% 

Galaxy Magnitude Error 

29% 

30% • 

31% 

Non-spinning Upper Limit (yT -1 !^ 1 ) 

8.7 x 1(T 3 

;2.2 X 10 3 

4.4 x 10~ 4 

Spinning Upper Limit (yr ” 1 L^ 1 ) 

! 2.7 x 10 -3 

5.3 x 10" 4 


TABLE I: Summary of results. The horizon distance is averaged over the time of the search. The cumulative luminosity combines the detection 
efficiercy with the galaxy catalog luminosity. Here, the value is the time-weighted average of the cumulative luminosity for each month. Many 
uncertainties are included in the calculation of the upper limit and they are summarized over all months. The effects of spin on BNS systems 
are negligible and not reported here. 




FIG. 4 The 90% rate upper limits as a function of mass. The first 
figure gives the upper limit on the rate of coalescence from BBH 
system as a function of the total mass of the system. The second 
figure gives the BHNS upper limit as a function of black hole mass, 
assuming a fixed neutron star mass of 1.35 M©. 


analyses, sky localization techniques are being developed to 
recons truct the direction of any gravitational- wave sources de- 
tected in the future. 
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Appendix A: Blind Injection Challenge 

During the blind injection challenge , simulated signals 
were inserted into the LIGO-only data without the search 
groups knowing the time or number of injections and their pa- 
rameters. Two blind injections occurred in the data described 
in this paper. The first simulated a burst of gravitational 
waves. The injected signal was the sum of two Gaussian mod- 
ulated sinusoids with linearly time- varying frequency. The 
root-square-sum amplitude h TSS for the signal was 1.0 x 10 -21 
at the Earth. The dominant component was at 58 Hz and the 
duration was about 12 ms. This injection was not a target of 
this analysis and was not identified as a significant candidate. 
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However, see Ref. [30] about the significance of this injection 
in the Collaboration’s burst search. 

The second blind injection was the simulated binary inspi- 
ral signal referred to in Section IV. The waveform simulated a 
binary system with masses 1.1 and 5.1 M©, with small spins 
0.19 and 0.06, respectively, in dimensionless units of the spin 
parameter a = ( cS ) / (Gm 2 ), at effective distance of 34,6 
Mpc for Hanford and 42.2 Mpc for Livingston. The candi- 
date identified in HI, H2, and LI has non-spinning templates 
with masses (1.0, 5.9), (1.0, 5.7), (1.1, 5.6) M©, and effective 
distances (43.6, 33.2, and 42.2) Mpc respectively. Given the 
parameters of the signal, the absence of a Virgo trigger in the 
coincidence does not cast any doubt on the validity of the can- 
didate. The loudest coincidence in the time slide background 
had a false alarm rate of 1 per 14 years. As this candidate 
was louder than that, 1 per 14 years is only a bound on its 
significance level. 

Since a candidate was identified coincident with this injec- 
tion, we conducted an extensive follow up study [23]. As part 
of the study, the SNR time series, x 2 time series and time- 
frequency spectrograms at the time of the candidate were in- 
spected. We also studied environmental influences on the de- 
tectors since this candidate was vetoed by a Category 3 data 
quality flag produced for high seismic noise at low frequen- 
cies at the Livingston Observatory. 

The 30-40 minute period of high seismic activity was due 


to earthquakes near Sumatra, which produced large ground 
motion at Livingston (but not Hanford) at frequencies between 
0.03 and 0. 1 Hz. There was a higher rate of accidental co- 
incidences in the time-slid data when using LI triggers pro- 
duced a few minutes after the time of the candidate. These 
accidental coincidences were coincident with peaks in seis- 
mic activity and excess power in the LI gravitational-wave 
channel. These accidental coincidences were H1L1 (double) 
coincidences, less significant than the triple coincidence can- 
didate seen at the time of the blind injection. No triple acci- 
dental coincidences were observed within that active seismic 
time. The time series of the gravitational -wave channel at the 
time of the candidate does not bear any resemblance to those 
at the times of the double coincidences correlated to seismic 
noise. 

The candidate passed all tests related to the pipeline and the 
statistical analysis. The presence of the seismic data quality 
flag in the time around this candidate does not substantially 
downgrade its significance. Had there not been a blind injec- 
tion at the time of this candidate, it would have been recog- 
nized as having an interesting level of statistical significance, 
and we would likely have pursued this candidate by dismiss- 
ing seismic activity as its cause. We are developing improved 
methods to better estimate the significance of such detection 
candidates. 
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